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Abstract. We show how the quasi- analytic method developed to solve linear elastic 
contacts to coated substrates (Perriot A. and Barthel E. J. Mat. Res., 2004, 19, 
600) may be extended to adhesive contacts. Substrate inhomogeneity lifts accidental 
degeneracies and highlights the general structure of the adhesive contact theory. We 
explicit the variation of the contact variables due to substrate inhomogeneity. The 
relation to other approaches based on Finite Element analysis is discussed. 



1. Introduction 

In elastic contact problems, it is known that the homogeneous substrate is a special case 
which leads to counterintuitive results. As an example, for adhesionless contacts, the 
relation between penetration S and contact radius a is independent from the mechanical 
parameters (Young's modulus E, Poisson ratio z/, reduced modulus E* = E/(l — u 2 )) 
and can be in all generality expressed as a function of the shape of the contacting bodies 
z(r) only. 

For instance, for a rigid cone of half included angle (3 in frictionless contact with 
an elastic substrate: 

vrtan/3 v ' 

Obviously, the reduced modulus is absent and one might conclude that the relation 
has a geometrical origin. To show that this is not the case, one needs only remove 
the assumption of homogeneous substrate. Then, for a coated substrate, for instance, 
the penetration will again depend upon the mechanical parameters of the system, in a 
highly non trivial manner pQ. Overlooking the more complex penetration dependence 
upon contact radius might actually result in inaccuracies in the determination of thin 
film mechanical properties by nanoindentation [21 (3J HJ [5] . 

Similarly, for adhesive contacts, a counterintuitive result is that the adhesion force 
or pull-out force F pu u out for a sphere of radius R and adhesion energy w is independent 
from the mechanical parameters. It is valued 

Fpullout = rnnuR (2) 

where various n values have been proposed such asl [6], 3/2 [7] or 2 [8]. 

| Corresponding author: etienne.barthel@saint-gobain.com 



Adhesive Contact to a Coated Elastic Substrate 



2 



However, changing the geometry from sphere to cone [H] or removing the assumption 
of a homogeneous substrate, for example with a coated substrate [TO], [11] , again results 
in reintroducing the dependence upon mechanical parameters. 

With the idea that the real structure of the elastic contact models only appears when 
considering a non homogeneous substrate, the present contribution aims at: 1) showing 
how the recently developed quasi analytical method for contacts to coated substrates [1] 
may be used to account for adhesion as previously used by Mary et al. [12] 2) exploring 
how the different concepts developed for homogeneous adhesive contacts may be applied 
to a non homogeneous system, to wit, a coated substrate. 

Section [2] reviews the various concepts which have been developed to establish the 
relation between macroscopic contact variables F, 5 and a and the physical phenomena 
of adhesion. Section [3] describes the surface stress transform method pioneered by 
Sneddon. Section @] explores the application of the method to the adhesive contact to a 
homogeneous substrate, while section [5] explores the application of the method to coated 
substrates. 

2. Adhesive elastic contact concepts and previous results on coated 
substrates 

2.1. Elastic adhesive contact - problem description 

Contact problems are characterized by mixed boundary conditions: the surface 
displacement is specified inside the contact zone and the normal stress at the surface 
outside. Note that this normal stress outside the contact zone is zero, except at the 
edge of the contact zone where more or less long range adhesive surface forces develop 
in the so-called interaction zone (Fig. [1]). 

Three concepts have been developed to handle these complex boundary conditions, 
provide a solution and establish relations between the macroscopic variables F, 5 and 
a. 

2.2. Three concepts to solve the adhesive elastic contact problem 

2.2.1. Energy minimization The first method disregards the contact edge details, 
assumes a given form for the surface stresses inside the contact zone and calculates the 
solution parameters by energy minimization, taking into account the adhesion energy 
wasa 7rwa 2 term. This macroscopic approach was implemented as early as 1934, when 
Derjaguin assumed a Hertz-like stress distribution and combined it with an adhesion 
energy contribution proportional to the contact area, which resulted [6] in a pull-out 
force given by n = 1 in Eq. [2j 

Based on elastomer adhesion experiments, a more realistic surface deformation built 
up from the addition of a flat punch displacement to the Hertz solution resulted in the 
JKR model [7] [jj where n = 3/2 in Eq. [21 This model, now generally accepted for soft 

§ Note that an early version of this model can be found in the appendix of G. Sperling's PhD thesis [13] . 
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solids, provides a more consistent description of the stresses inside the contact zone than 
the original 1934 result. 

More generally, a formal method to minimize the total energy with regards to 
the contact radius a naturally leads to introducing the energy release rate, a concept 
borrowed from fracture mechanics : if S is the stored elastic energy, the energy release 
rate G can be defined in an axisymmetric geometry as 



dS 

2-naG(a) = — — 
da 



(3) 

5 

Such an approach was used for instance as an alternative derivation of the JKR model 
in rE 



2.2.2. Stress intensity factor The second method is to focus on the crack-like stress 
singularity at the contact edge. It is described in terms of the stress intensity factor. 
The relation between the energy release rate and the stress intensity factor was initially 
proposed by Irwin [15] . In the field of the adhesive contact, this second method was 
mainly propounded by Maugis [91 HH [T7] . It is particularly useful for direct extensions 
of the Linear Elastic Fracture Mechanics where a description in terms of stresses at the 
crack tip is necessary: for instance, it may be used to include crack tip viscoelasticity 

2.2.3. Exact solution The third method is actually more general: it takes into account 
the details of the attractive stress distribution outside the contact zone. This idea 
was somehow initiated by the DMT model [8] (n = 2 in Eq. [2]) but was brought to 
a new dimension by the introduction of a cohesive zone model by Maugis [19], which 
resulted in a clear exposition of the transition between the two limiting models JKR 
and DMT. It also led to the development of the fully viscoelastic adhesive contact 
models [201 [21] . In these models, stress relaxation inside the contact zone - as in the 
adhesionless viscoelastic contact by Ting [22J - and creep inside the interaction zone - 
as in the viscoelastic crack models [TBI El HE] - are taken into account simultaneously. 

Just as the viscoelastic case extends the adhesive contact model from instantaneous 
response to time dependent constitutive relation, the coated substrate generalizes the 
model to an inhomogeneous half space. Our interest here will however be limited to 
the small cohesive zone sizes, located on the JKR side of the transition, which connect 
this approach to the previous two. For an extension to the JKR-DMT transition with 
a coated substrate, when the interaction zone grows to extensions comparable or larger 
than the film thickness, in the spirit of the Maugis model, see [23] . 

where it stands as a curiosity, not as the "Theorie der Haftung" developped in the body of the text, 
which mainly deals with surface roughness and plastic deformation. 
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2. 3. Review of previous results on the adhesive contact to coated substrates 

On the theoretical side, the literature on the adhesive contact to coated substrates 
is scanty, although, in practice, thin soft adhesive layers are often used, as in many 
applications of pressure sensitive adhesives and functional or protective organic coatings. 
Actually only two series of papers by Shull and coworkers [2H [25] on one side and 
Johnson and Sridhar [TQl EH] on the other side have been published. 

Shull and coworkers calculate the energy release rate from the contact compliance 
(cf. I5.2.2p . To complete the calculations, the compliance of coated sustrates are 
computed by Finite Element (FE) methods and phenomenological fits or numerical 
differentiation may be carried out. 

Johnson and Sridhar's approach goes along the line of the stress intensity factor 
(section 12.2.2ft approach: still using FE, they directly calculate the energy release rate 
G by a stiffness matrix variation method [26] and extract stress intensity factors from 
G. From this stress intensity factor, they build the adhesive contact relations. 

Both methods converge on their ultimate use of Finite Element calculations to 
handle the complexities of the coated substrate response. Appearances to the contrary, 
these two series of papers seem actually to rely on the same method and an in depth 
investigation of the Parks method in the context of the frictionless contact would 
probably connect the stiffness matrix derivative used by Parks to the contact stiffness 
derivative method as developped by Shull. 

In the present paper, we show that the quasi-analytical method introduced 
previously pQ may be used to efficiently compute all the necessary contact variables 
and in particular provide direct numerical evaluations for the stress intensity factor for 
coated substrates. It is also useful to expose the equivalence of the various concepts for 
the adhesive contact in the wider context of non-homogeneous substrates. 

3. Sneddon and the surface stress transform 

Adequate transforms facilitate the investigation of linear elastic axisymmetric contacts 
formulated as integral equations [27] . It is ideally suited to the adhesive contacts because 
of the relevant boundary conditions including the interaction zone. 

Following our previous studies [20], [21] [28] [29], we use g(r), which is a transform of 
the normal stress distribution at the surface a z defined as 



and 9(r) which is a transform of the normal surface displacement u z defined by 



These relations are suited to the adhesive contact problem because g(r) is expressed 
as a function of normal surface stresses at radii values larger than r and 9(r) as a 
function of surface normal displacement at radii values smaller than r, in agreement 




(4) 




(5) 
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with the adhesive contact boundary conditions (Fig. [Q. In particular g is zero outside 
the interaction zone. 

Similarly, inside the contact zone (r < a), u z (r) = 5 — z(r) where z(r) is the shape 
of the indenter. Integration by parts transforms Eq. [5] into 

9(r) = 5-5 (r). (6) 

where 5 is the penetration and 5q depends only upon the shape of the indenter z(r) 
through 

. . s f r z'(s)ds 
Jo sjr 1 — s z 

The full power of these transforms appears when it is recognized [20l |2T] that in 
the linear elastic case and for a homogeneous substrate, mechanical equilibrium results 
in 

9(r) = ^0(r) (8) 
This provides the direct solution to the adhesive contact problem which we detail 

now. 



4. Adhesive contacts on homogeneous substrates 

Before tackling the problem of the adhesion to coated substrates we review our 
present understanding of elastic contacts to homogeneous substrates using the transform 
method. 



4-1. Non adhesive contact of smooth indenters 

In the absence of adhesion, Eq. H] shows that g(a) = 0. If the indenter shape z(r) is 
smooth at the contact edge, then Eqs. [6] and [8] imply 8(a) = and 

S = 5 (a) (9) 

Thus the indenter penetration in an adhesionless contact is the function 60(a), 
which depends only upon the geometry of the indenter (Eq. [7]) and not on the elastic 
parameters. This is the counterintuitive result stated in section [TJ 

4-2. Non adhesive contact of flat punches 

A special case however is the flat punch for which £o( r ) = for r < a and g(r) = 
for r > a. The penetration is independent of the contact radius as expected since the 
contact radius is constant, equal to the punch radius. 

Then 9(a~) = 5 so that g(a + ) = and g(a~) 7^ 0. The exact meaning of this 
discontinuity in the g transform at a will be discussed below (section [4.4.21) . 

Note also that the force can be obtained through the simple but general 
expression [29~] 

r+00 

F = 4 / g(s)ds (10) 
Jo 



Adhesive Contact to a Coated Elastic Substrate 



6 



Specifically here, 

F (a) =4 P g{s)ds (11) 
Jo 

For the flat punch the force is then directly calculated from Eqs. [6] [8] and [10] as 

F (s) = S(a)8 fp (12) 

with 

S(a) = 2aE\ (13) 
4-3. Adhesive contacts 

In the adhesive case, however, although the indenter shape is smooth at the contact 
edges, the condition g(a) = is relaxed due to the adhesive interactions. The core 
of the JKR approximation is to neglect the details of the stresses and deformations at 
the contact edge and encapsulate the adhesive contribution in an additional flat punch 
displacement. This amounts to a description where the attractive interaction stresses 
drop abruptly to zero outside the contact zone, leading to a stress singularity akin to 
those met in fracture mechanics. 
Then Eqs. [9] and [11] become 

6 =6 + 6 fp (14) 

F — F + S(a)5 fp (15) 

The question is to provide a relation between this additional flat punch displacement 
8f p and the adhesion energy w for a given contact radius a. 

4-4- Various concepts for the derivation of adhesive contact results 

For future reference, we now apply the surface stress transform solution to the three 
methods developed to handle adhesive contacts (section |2~2]) . 

4-4- 1- Energy Release Rate A very compact derivation is obtained when it is recognized 
that for a homogeneous substrate g et 9 are independent of a inside the contact zone. 
In addition the total mechanical energy £ as a function of the normal surface stress 
distribution is [29] 

A r+oo r+oo 

£ = — ds g 2 (s) = E* ds 9 2 (s) (16) 
E* Jo Jo 

Then for constant penetration, with g(a) = outside the contact zone, 

2naG(a) = ^4r^ = E*6(a) 2 (17) 
E* 

Equating G and the adhesion energy w results in 

« = ^ (18) 
wE*a V ; 

The energy release rate G(a) can be expressed as a function of the local stress 
distribution at the crack tip because of Eq. [Hand g(r) = for r > a. 
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4-4-2. Stress intensity factor This expression together with the plain strain Irwin 
relation 0, [15] between stress intensity factor and energy release ratqjjj 

G = K I 2 /{2E k ) (19) 

suggests that 2g(a)/ \fjra assumes the role of a stress intensity factor. 

A direct derivation is obtained [9] if we use the following expression for the stress 
intensity factor at the contact edge 

K l = lim (2ne) 1/2 a(a - e) (20) 
Eq. H] can be inverted [29] providing 



( \ 2 
o- z {s) = - 

7T 



g(a)Q(a-s) + f +QO df _ 9'(t) 



(a^-s 2 ) 1 / 2 ' J s (t 2 - s 2 ) 1 / 3 ^ 
where it has been assumed that g is smooth except for a discontinuity at a. We use the 
notation g(a~) = g(a) and assume g(a + ) = 0. Then 

K, = ^ (22) 
\/7ra 



which confirms that g(a) has the form and meaning of a stress intensity factor. 

As a stress intensity factor, g(a) is a measure of stress singularity at the contact 
edge. An ancillary property of the transform defined by Eq. H] is that it regularizes 
the singular crack like stress distributions and Eq. [5] establishes a relation between the 
additional flat punch penetration and the stress singularity through 

0(a) = 6 fp = 2 -^ (23) 

4-4-3- Self- consistent method The self-consistent description of the interaction zone 
[T9| l28l [30] can be explored at the limit of negligible interaction zone extension, i.e. in 
the JKR limit. 

One possible starting point of the self-consistent method [281 EH] is to calculate the 
adhesion energy by 

r+°° , dh , . , 

w = - a(s)—ds 24 

Ja ds 

Following [21] (Eq 15) the gap between the surfaces h is split into the contributions of 
the contact stresses and the interaction stresses 

h(r) = h Hertz (r, a) + h int (r, a) (25) 

If the interaction range is small (which results in c — a << a), the radial extension of the 
interaction will also be small, g' is peaked around a, and the dominant term in Eq. [2H 
will be the du z ^ nt /ds term. Thus |29j, 

As a result, one again recovers Eq. [THJ Combined with Eq. [22], this method may be 
viewed as a direct derivation of the Irwin relation |15j . 

|| The 1/2 factor comes from the rigidity of the punch. 
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4 ■ 5. Normalized form 

The adhesive contact equations are obtained by combining Eqs. [T4lfl5l and Eq. [T8j For 
later reference, we specialize the results to the spherical indenter and use the Maugis 
normalization [19] 



P = ^B (27) 
ttwR 

D = 6 — m (28) 



( 7T 2 U> : 
E* 



(29) 



1/3 



One obtains 



P = -a 3 - 2 v / 2a 3/2 (30) 
3 

D = a 2 - V2a 1/2 (31) 

These are the usual JKR equations which will be generalized in the next section. Their 
structure directly reflects the general structure of Eqs. [T4TfT5l with the identifications 
5f p = —\[2a i l 2 and Sf p = 2a. 



5. Coated Substrates — Thin Films 

5.1. Description of the contact 

If the substrate is not homogeneous the simple diagonal equilibrium relation Eq. [8] 
between g and 6 is lost. However, keeping the same transforms, we have shown pQ that 
a useful relation subsists 

d(r; a, t, [E]) = 2g ^ a ^ ^ + I f + °° g ( s; a , t, [E])K(r, s; t, [E})ds (32) 

Ei IT Jo 

where K expresses the elastic response of the coated substrate [I] , t stands for the coating 
thickness and [E] for the four mechanical parameters (film and substrate Young's moduli 
and Poisson ratios). Since g is zero outside the contact zone, the upper boundary will 
typically be the contact radius a or the interaction zone radius c if c 7^ a. The finite 
integration interval facilitates the numerical inversion [1] of Eq. [321 Note that direct 
analytical inversion is impossible because of the complex expression for K. 

We now split the stress distribution and the displacement field into their non 
adhesive 5h and flat punch 5f p components. 

For the non adhesive contact of an indenter of regular shape, the penetration is 
given by 

x I + \m\ x / \ 2 9H(r;a,t, [E]) 

5 H {a, t, [E]) - 8 [r) = — (33) 

+ - f 9h(s; a, t, [E))K{r, s; t, [E])ds (34) 

IT JO 
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where 8 (r) reflects the indenter shape and <7#(a; a, t, [E]) = in the absence of adhesion, 
ensuring a unique solution. Eq. [3H generalizes Eq. [15] and introduces the non-trivial 
relation mentioned in section [U 

For the contact of the flat punch the problem is linear with 5f p and 

= 2g fp {r;a,t, [E}) + 2 r« ^ ^ ^ ^ ^ 

Hi i 7T JO 

with gf p (r; a,t, [E]) oc 5f p and g(a; a,t, [E]) ^ 0. Eq. [35] generalizes Eq. [23] into a non 
trivial proportionality relation between g(a) and the flat punch penetration 5f p . 

Then the penetration for the adhesive contact of an indenter of arbitrary shape is 

8{a,t,[E]) = 6 H {a,t,[E]) + 8 fp (36) 

With the following definition of the contact stiffness 

S(a,t,[E})= Ffpia : t,m (37) 
°fp 

the force is 

F(a,t, [E]) = F H (a,t, [E]) + S(a,t, [E])5 fp (38) 

Equations [36] and EE] generalize Eqs. [15] and [HI Here again, the question is to 
provide a relation between the additional flat punch displacement 5f p and the adhesive 
interaction. 



5.2. Derivation 

For coated substrates, the resolution benefits from more general expressions derived for 
viscoelastic adhesive contacts [31]. This is not fortuitous but results from a similar 
breakdown of the simplified relations derived through Eq. [8] when spatial heterogeneity 
or time dependence is introduced. 



5.2.1. Energy Release Rate The total energy is [31J 

r+oo 

S = 2 ds g{s)6{s) (39) 



Then 

2naG(a) = 2 (g(a)9(a) + £ ds ^j^(s)) (40) 

because G(a) is calculated at constant displacement so that inside the contact zone, 
the surface displacement - and therefore 6(r) - are unaffected by the additional stress 
distribution g(a): there is only a stress rearrangement inside the contact zone which 
does affect g(r), r < a for a coated substrate: indeed the local relation Eq. [8] between g 
and 9 breaks down in this case. 

Following Mary et al. [12] , one can show that this non local contribution is actually 
canceled by the non local contribution included in the g(a)6(a) term. For s < a 

n d9{s) -g(a)K(a,s)+ f dr ^K(r,s) (41) 



Ei* da Jo da 
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and multiplying by g(s) and integrating between and a, one obtains 

The region affected by the crack tip stresses extend over a distance commensurate 
with the interaction zone size, i.e. it is smaller than the film thickness so that the energy 
release rate is controlled by the film compliance. 

5.2.2. Compliance Method The compliance formulation is at the core of the method 
used by Shull and coworkers but the derivation given in some of their earlier papers 
was obscured by unnecessary assumptions. This formulation emerges readily when £ in 
Eq. [39] is expressed as a function of strain instead of stress^- 

We split the total stress and strain fields in their non adhesive (gn and 9h) 
and flat punch components (gf p and 9f p ). This allows for easy integration because 
9f p (r;a,t, [E]) = 5f p is independent of r. Using Betti's theorem to calculate the cross 
terms 

r+oo 

2 dsg H {s;a,t,[E])9 fp {s;a,t,[E]) (43) 

J 

r+oo 

= 2 dsg fp {s;a,t,[E})9 H {s;a,t,[E}) (44) 

J 

= l ~8 fp P H {a,t,[E\) (45) 
and similarly using Eq. I3T1 

r+oo \ 

2 J o dsg fp {s;a,t, [E])9 fp {s;a,t, [E]) = -S(a,t, [E})8% (46) 

to express the flat punch elastic energy, one obtains the total energy 

£(a,5 fp ) = U H (a,t, [E]) + ~S(a,t, [E\)8% + P H (a,t, [E])5 fp (47) 

where 5jp is negative. 

A graphic illustration of this result is given in Fig. [21 

The energy release rate G is the differential of the total energy S(a, 5f p ) with respect 
to contact area (Eq. [3]) at constant total penetration S where 

5 = 8 H (a) + 5 fp (48) 



dS 
da 



d£ ^ d£ d5f p 
da d5f p da 



Now 



dU H = P H (a,t,[E})dS H (50) 
dP H = S{a,t, [E])d8„ (51) 

Such a strain energy release rate has also been introduced previously in viscoelastic crack problems 
by Schapery [32] and Greenwood and Johnson [18] , 
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so that finally all terms cancel^] except the differential of the flat punch elastic energy 
and 

2naG(a) = ±8%^- (52) 
This differentiation process is also illustrated on Fig. [2j 

5.2.3. Stress Intensity Factor The derivation of the stress intensity factor expression in 
section 14.4.21 is indendent of the material properties and only results from the definition 
of the transform g. It is therefore unchanged in the more general case and the pending 
problem is actually the relation between g(a) and the energy release rate G. 

5.2.4- Self- consistency In Eq. |25"1 

, . . 2 r 0(s) - 0(a) , , , 

h mt (r, a, t,[E]) = - K L^ } ds (53) 

If at the crack tip the deformation predominantly results from the adhesive interactions 
- and this is the essence of the JKR limit - then Eq. [2H becomes [20] 

2 r c 

w = / 9'(r)g(r)dr (54) 

na J a 

From Eq. [32] we obtain 

, (r ) W + Y" s(s) f^£) (55) 
ti\ 7r Jo ar 

The second term is well behaved when c — > a since K is the elastic response of the 

coated substrate: the contact edge singularity results from the boundary conditions, 

not the response function. Thus it is g' which is peaked around a [28] so that in the end 

for c < r < a 

6>{r) * ^ (56) 

and again 

2fl(a) 2 r , 7 , 

irE*a 

even for a non-homogeneous substrate. 

The local response is indeed dominated by the local stress distribution and the 
film compliance: the same relation Eq. [T7| between adhesion energy and stress intensity 
factor applies. 

+ Due to the counterintuitive results obtained for the homogeneous substrate, one may develop doubts 
about the identity of the stiffness defined by Eq. [51] and the stiffness defined by Eq. [37l That this 
identity does hold for a non homogeneous substrate is shown in the appendix (section [8] ft is due to 
the absence of adhesion: the normal stress at the edge of the contact is zero so that contact radius 
variation does not result in force variation, to first order. 
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5.3. Relation between g (a) and S 

Comparing Eqs. [52] and [57] one infers 
Ag{af _ 1 2 dS 

~ 2^Ta ^ 
This is a remarkable result which one should in principle be able to derive from Eq. 

Note however that the present quasi-analytical method [I] simultaneously provides 
the flat punch stiffness and and contact edge stress singularity g(a) (or the stiffness 
derivative through Eq. [35]) by the simple resolution of a linear system. 

5.4- Normalized Solution 

The procedure to compute a full force curve for the adhesive contact to coated substrates 
is therefore: for a given value of the contact radius a 

(i) compute the adhesionless penetration and force for the given indenter shape 

(ii) compute the force and stress intensity factor (or g(a)) for the flat punch for a unit 
value of the penetration 

(iii) compute actual value of g(a) from Eq. [57] and rescale flat punch force and 
penetration 

(iv) compute the solution from Eqs. [36] and EE] 

This provides the solution under the form of two relations between the local 
deformation (stress intensity factor) and the local response on the one hand and the 
remote loading and the macroscopic response (contact stiffness) on the other hand. 

In normalised form, with 

r 



one haa3 



t 



(59) 



Fo, s =^|-n.(a, *,[£]) (60) 

So,s = ^A s (~a,t,[E}) (61) 

S = 2aE 1 *£ eq (d, t, [E]) (62) 

g(r-a,t,[E]) = 5 -^T(r-d,t,[E]) (63) 

where the normalized variables can be numerically calculated by the simple 
algorithm presented in [1]. 



For a cone: 



Fo ' c =4W nc(M ' [i?]) 
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All variables equal 1 for the homogeneous substrate except Tl s = 8/3 for the spher 
r(f; a, t, [E]) is the surface stress transform normalized to penetration 5. In 
particular, T(d; a,t, [E]), denoted below for brevity, is the contact edge singularity 
g(a; a, t, [E]) incurred for a coated system - normalized to a homogeneous material with 
the film elastic properties - at identical 5f p value. The variable which is positive 
since both 5f p and g(a) are negative in Eq. [631 is a function of a, and depends upon the 
mechanical parameters of the system. From Eq. [52] the following identity holds: 



r(i) = ./-L^ = J«&J (64) 

V ; y 2ES da V da K ' 

Then for the sphere, keeping the Maugis normalization, one introduces the film 
thickness normalized to the typical adhesive contact radius 

v = Wf" (65) 

and with a = rjd one obtains 

n s=W 3 ^-2V2W 3/2 ^ (66) 

D s = ( V ~a) 2 A sfi - V2 (rjdf 2 ^- (67) 

These equations generalize Eqs. I30T13T1 The homogeneous substrate force and 
penetration terms are corrected by the coated system factors n Sj0 and A^o, the 
homogeneous contact stiffness by £ cq and, for identical stress intensity factor, the 
homogeneous penetration is corrected by 1/T(1). 



6. Examples of Numerical Results 

6.1. Reduced variables 

For given mechanical parameters Ex*/Eq*, vq and v\ one may calculate the four variables 
n Si o,A Si0 ,^eq and T(l) as a function of a = a/t. Typical results are illustrated in Fig. [3] 
for Ei* /Eq* = 10 and Fig. H]for E 1 */E * = 0.1. The results for these reduced variables 
compare well with the FE calculations by Sridhar and Johnson [Til . 

The thin film contact problem exhibits a transition between film dominated to 
substrate dominated behaviour. In each limit cases, the contact behaves like a contact 
to a homogeneous system. The non trivial behaviour is apparent in the transition which, 
roughly speaking, occurs for a ~ 1, but is shifted to higher values for compliant films 
(Fig. [3]) and to lower values for stiff films (Fig. HJ) [TJ. 

£ eq exhibits a transition between film to substrate reduced modulus as the contact 
radius increases is consistent with numerous previous works [33]. The behaviour of A S; o, 
which deviates from 1 in the midst of the transition |1] |5] has been much less studied, as 

jj For a cone Il c = 2 
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mentioned previously (section [T]) . Similarly for the flat punch we note that T(l) tracks 
S eq at small contact radius values, up to about the middle of the transition. This does 
result from Eq. EH as a linear expansion of £ eq shows (note that £ eq (0) = 1). For larger 
contact radius values, higher order terms come into play and T(l) starts to deviate from 

'-'eq- 

6.2. Adhesive contact solutions 

The reduced variables II s (a, t, [E]), A s (a,t, [E]), £ eq (a,t, [E]) and T(l) = T(a;a,t, [E]) 
provide the normalized solution (Eqs. 1661 and 1671) as a function of 77 for a set of mechanical 
parameters. It is worth emphasizing that - provided the reduced variables have been 
calculated on a wide enough range of values of a = a/t - the solution for any values of 
contact radius, film thickness and adhesion energy can be calculated from the reduced 
variables through the elementary arithmetics of the normalization Eqs. [66] and [671 This 
should provide for an easy algorithm to fit experimental data for adhesive contact on 
thin films. 

Results for various values of r\ are displayed on Fig. [5] and [6] for E\ / E * = 10 and 
Fig. [7J for E^/Eq* = 0.1 with v = v x = 0.25. 

For a coated substrate in the presence of adhesion, the system undergoes two 
transitions. The first one is is the transition between coating and substrate response p], 
when a ~ 1 as described in the previous section. The second one is the transition from 
adhesion dominated to non adhesive contact as the load and therefore the contact radius 
increase [31]. This transition occurs when a ~ 1. The deviation from the JKR results 
occur when the two transitions are simultaneous, i.e. when 77 ~ 1. 

Explicitly, for large values of rj (i.e. thick films), the adhesive stage of the contact 
will occur for small a and the contact is film dominated and behaves as a homogeneous 
material with the mechanical parameters of the film. 

For small values of rj (i.e. thick films), the adhesive stage occurs for large a and the 
contact is substrate dominated: the contact behaves as a homogeneous material with the 
mechanical parameters of the substrate. However, at the edge of the contact, this contact 
is not equivalent to an adhesive contact to a homogeneous material with substrate 
properties. Indeed, the stress intensity factor is dominated by the film properties 
(Eq.ET]). 

As a result, moderate deviations on the pull-out force are evidenced. For soft layers, 
the adhesion force is enhanced [II] . For rigid layers there is a small reduction of the 
adhesion forc iffl . 

ff This is not, however, the reason why on soft materials a rigid layer may drastically reduce the 
adhesion, which is often due to the suppression of non elastic additional dissipation phenomena. Note 
that here the adhesive force reduction is accompanied by an increase of the stress intensity factor. 
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7. Conclusion 



The surface stress transform is adequate to handle the adhesive contact to coated 
substrates. It provides a numerically simple method to compute the four quantities 
necessary for the actual description of such adhesive contacts. For a given layer and 
substrate mechanical properties, fits to data with free adhesion energies and coating 
thicknesses may be performed easily. It also allows a consistent description of the crack 
tip which connects the various concepts developed in adhesive contact problems. 

The method also allows the inclusion of finite indenter stiffness, a question which 
arises in practice and which will be studied in more details in a subsequent paper. 
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8. Appendix: contact stiffness and flat punch stiffness 

The identity of the contact stiffnesses defined by Eq. [37] and by Eq. [51] may look 
questionable: is the contact stiffness for an adhesionless curved indenter identical with 
the flat punch stiffness for an identical contact radius, even for an inhomogeneous 
substrate ? 

Since the curved indenter is adhesionless, g(a) = 0. Then from Eq. [10] we have 

dP r a dg(r) . , na . 
do Jo do 

From the differentiation of Eq. [32] with 5 together with Eq. [5] 

1 = - r^P-K(r,s;t,[E])dr (69) 

7T JO do 

This is the flat punch equilibrium equation for unit penetration. Then the force 

— (TO, 

is the flat punch force for unit penetration, i.e. the flat punch stiffness S(a). Therefore 
for an adhesionless curved indenter 

f = s(a) (71) 

Note however that in general this stiffness is not the hertzian 3PR/4a 3 . 
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Captions 

Fig. 1: schematics of an adhesive contact: the contact zone (radius a) is surrounded 
by the interaction zone. The gap between surfaces outside of the contact zone is h(r). 
The normal stress distribution turns from compressive at the center of the contact to 
tensile on both sides of the contact edge. The transform 9(r±) is calculated from the 
normal surface displacement at r < r x . The transform g(r 2 ) is calculated from the 
normal surface stress for r > r 2 . 

Fig. 2: Adhesive contact configuration obtained from the adhesionless contact by a 
flat punch displacement 5f p along the tangent to the adhesionless force F-penetration S 
curve. Integration gives the stored elastic energy. The energy release rate is obtained 
as the variation with contact radius a of the stored elastic energy at constant total 
penetration 5 (shaded 

Fig. 3: Normalized force IT S (left) and penetration A s (right) for an adhesionless 
sphere as a function of contact radius normalized to film thickness a = a/t. The film to 
substrate modulus ratio is 0.1 and Poisson ratio 0.25 for both materials. The contact 
stiffness E eq (left) and contact edge stress intensity factor T(l) (left) are calculated from 
the flat punch boundary conditions. For a ~ 1 the transition from film dominated to 
substrate dominated contact occurs. 

Fig. 4: Same plot as Fig. 3 for a film to substrate modulus ratio equal to 10. The 
transition occurs earlier due to film rigidity. 

Fig. 5: Normalized force versus normalized contact radius as a function of r), the 
ratio of the film thickness to the typical zero load adhesive contact radius. When 
7] ~ 1, adhesive effects are significant in the transition from film to substrate dominated 
regimes: deviations from the JKR model are observed and increased adhesion forces are 
observed. 

Fig. 6: Normalized force versus normalized penetration as a function of r\. Identical 
parameters as in Fig. 5. 

Fig. 7: Similar plot as in Fig. 5 with a film to substrate modulus ratio of 10. In 
the intermediate regime r\ ~ 1, adhesion forces are reduced. 
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